Model of the Influence of an External Magnetic Field on the Gain of Terahertz 
Radiation from Semiconductor Superlattices 
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We theoretically analyze the influence of magnetic fleld on small-signal absorption and gain in a 
superlattice. We predict a very large and tunable THz gain due to nonhnear cyclotron oscillations 
in crossed electric and magnetic fields. In contrast to Bloch gain, here the superlattice is in an 
electrically stable state. We also find that THz Bloch gain can be significantly enhanced with a 
perpendicular magnetic fleld. If the magnetic fleld is tilted with respect to the superlattice axis, the 
usually unstable Bloch gain profile becomes stable in the vicinity of Stark-cyclotron resonances. 
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Semiconductor superlattice (SL) is a model system for 
a wealth of fundamental phenomena resulting from the 
wave-nature of charge carriers [li [3, S| ■ Interesting ex- 
amples of such phenomena are Bloch oscillations in dc 
biased SLs [3| , Shapiro- like steps Q and parametric res- 
onance in time-dependent electric fields, coherent Hall 
effect in crossed electric and magnetic fields and Stark- 
cyclotron resonances in tilted magnetic field The 
semiclassical theory predicts that electrons performing 
Bloch oscillations in the presence of weak dissipation can 
potentially provide THz Bloch gain ^9!|. The Bloch gain 
profile, which is shaped as a familiar dispersion curve, is 
not limited to SLs. It was recently predicted [l^ and ob- 
served TT] in intersubband transitions of quantum cas- 
cade lasers, but dispersive gain profiles have been also 
found in the microwave responses of Josephson junctions 
and Thim amplifiers ^is'l . Here we consider the Bloch 
gain in its traditional meaning as an effect occurring due 
to Bloch oscillations in a single energy band. In SLs, the 
realization of Bloch oscillator is a long-standing prob- 
lem due to the instability of a homogeneous electric field 
in conditions of negative differential conductivity (NDC) 
14|. This electric instability results in a formation of 
electric domains in long SLs, which are destructive for 
the THz gain. 

A considerable amount of theoretical and experimental 
activities have been devoted to the investigations of tran- 
sient oscillations and voltage-current (VI) characteristics 
in the presence of magnetic field [1, 0, 0, 0, llfi, 17, 18], 
but the response to a time-dependent field is mostly un- 
explored. In this letter, we focus on the influence of 
the magnetic field on the small-signal absorption and 
gain. We find that absorption and gain profiles in crossed 
electric and magnetic fields have different characteristic 
shapes in the Bloch-like and cyclotron- like regimes of bal- 
listic motion. We predict a very large and tunable THz 
gain due to nonlinear cyclotron oscillations in a single 
energy band. The cyclotron gain takes place in con- 
ditions of positive differential conductivity (PDC) and 
therefore it is stable against space-charge fluctuations. 



We also demonstrate that Bloch gain can be significantly 
enhanced by the perpendicular magnetic field. Moreover, 
the usually unstable Bloch gain profile becomes stable in 
the vicinity of Stark-cyclotron resonances if the magnetic 
field is tilted with respect to the SL axis. 

We consider a SL under an action of a static magnetic 
field B in arbitrary direction and an electric field E{t) — 
Edc+Ei^ cos Lut in a SL direction, which is chosen to be the 
x-direction. Here iJ^j cos cut is a weak probe field with a 
frequency co fixed by an external circuit (external cavity) . 
We consider the electron transport in a single miniband 
with the standard tight-binding dispersion relation 

e(k) = -{A/2)cosk^d+h^{kl + kl)/2m, (1) 

where e(k) is the electron energy, k is its quasimomen- 
tum, A is the miniband width, d is the SL period and m 
is the effective mass in y- and z-directions. We use the 
semiclassical approach based on the Boltzmann equation 
in the relaxation time r approximation 
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where Vi{ki) = h~^de(k)/dki (i = x,y,z) are the velocity 
components and feq (k) is the Fermi distribution 0, 0] . 
By solving Eq. ^ , we find the stationary time-dependent 
current after transient 
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d^fc /,g(k) / ds e-(*-^)/^v(k*), (3) 



where the prefactor takes into account the density of 
states [11, kx is integrated over the Brillouin zone and 
the integration limits for ky and fc^ are zLco. Here k* is 
a ballistic trajectory governed by the equations 
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eE(t) + ev(k*) x B 
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The second equation in (jj]) means that the quasimomen- 
tum at i = s is k. We derived the solution of Boltz- 
mann equation ([3]) by using both a generalization of the 
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technique based on time-evolution operator [18| and the 
method of characteristics [l^. The approach based on 
this solution is most powerful in the limits of low carrier 
density N and low temperature. In this case 



(5) 



and the current ([3]) is determined by the ballistic trajec- 
tories starting at k = 0. 

The time-dependent current ([3]) can be used to calcu- 
late the real part of the dynamical conductivity ariLo) = 
Re[cr(cL;)], which determines the gain [ar < 0) and ab- 
sorption (CTr > 0). In the absence of magnetic field ar is 
defined by the Tucker formula 0, [lO] 

jdc{eEdcd+hu}) - jdc{eEdcd- hhj) , 
aAu^) = — ed, (6) 

where jdci^Edcd) is the Esaki-Tsu characteristic fl*], 

2eEdcdT/h 



jdc{eEdcd) = jp 



1 -I- [eEdcdr/hf 



(7) 



n = 
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and jp is the peak current corresponding to the critical 
field Ecr = h/edr. The Drudc conductivity of the SL 
is (To — 2jp/Ecr- If Edc > Ecr, Eq- & describes the 
dispersive Bloch gain profile [blue line in Fig. [3] (b)] with 
a crossover from gain to loss at the resonance uj k lob, 
where lob — eEdcd/h is the Bloch frequency. As directly 
follows from Eqs. ^ and ([7]), the upper limit for this 
Bloch gain is min[(Tr(ti-')] = — ctq/S. 

We now turn to the consideration of electron dynam- 
ics in crossed electric and magnetic fields B = (0, 0, Bz). 
The ballistic electron dynamics under static fields is de- 
termined by the pendulum equation, which follows from 
Eq. Q 0, Il5i . If the electrons start at the bottom of the 
miniband, the frequency of their nonlinear oscillations is 

■mjJc/1K{LJB I'i.LJc), wb < 2aJc 
ttujbI'2.K{2uJc/lob)., ujb > '2lOc, 

where ujc — gB^ j ^nrixm and rrix — 2h^ j /S.d^ arc the cy- 
clotron frequency and the effective mass at the bottom 
of the miniband, respectively. Here K{k) is the com- 
plete elliptic integral of the first kind as a function of 
elliptic modulus. There exist two distinct regimes of mo- 
tion: cyclotron-like oscillations for large magnetic field 
ujb < (oscillations of pendulum) and Bloch-like os- 
cillations for dominating electric field lub > 2u!c (rota- 
tions of pendulum). As can be seen from Fig. [1] (a), 
the frequency f2 is tunable by variation of the magnetic 
and electric fields. If ujb ^ i^o we get from Eq. ^ 
that Q = UJB- In the opposite limit ^ we have 
= uJc- Close to the separatrix lub = 2uJc oscillations 
are strongly nonlinear and anharmonic. All harmonics 
of f2 are present in the Bloch-like regime whereas only 
odd harmonics exist in the cyclotron-like regime. Such 
kind of transient oscillations have been directly observed 
in the experiments [7[. 
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FIG. 1: (color online), (a) Frequency of the nonlinear oscil- 
lations n [Eq. ((Sjl] as a function of lob for uJcT 2,4. The 
dashed lines separate the cyclotron- and Bloch-like regimes 
of motion, (b) VI characteristics for uicT =: 0,2,4. (c)-(f) 
Absorption and gain profiles for ujcT — 4 and different values 
of electric field ujet (c), 6 (d), 7.5 (e) and 8.5 (f). 



In the presence of scattering these oscillations deter- 
mine the dc current density at low temperatures, as can 
be directly seen from Eqs. dSl)-® [IBl- The resulting VI 
characteristics for different values of the magnetic field 
are shown in Fig. [1] (b). The transition from oscillatory 
to rotational motion at the separatrix ujb = 2lUc mani- 
fests itself as a very abrupt and strong suppression of the 
dc current 3, 3 12 1- In the cyclotron-like regime, there 
are no Bragg reflections and therefore it corresponds the 
PDC part of the VI characteristic. On the other hand, 
the Bragg reflections in the Bloch-like regime result in 
NDC at the operation point. We find that the response 
of the miniband electrons to a weak ac field also shows 
clear signatures of the different types of ballistic motion. 
Figs.[l](c)-(f) show the absorption and gain profiles, cal- 
culated using Eqs. ([3])-([5]), for different values of ojc and 
UJB- The field strengths are chosen in such a way that the 
electrons perform several cycles of oscillations between 
the scattering events J7r > 1. We see that the gain and 
absorption profiles in the different regimes of oscillations 
have their own characteristic shapes. If Wc ^ ^b, elec- 
trons are restricted to the bottom of the miniband, and 
we obtain a familiar Lorentzian absorption profile of a 
harmonic oscillator, which is centered at w « [Fig. [1] 
(c)]. In this case ar is always positive and an amplifica- 
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FIG. 2: (color online), (a) Cyclotron gain profiles for wbt — 
3.6, LOcT — 2 and lubt — 8.5, ojcT = 5. (b) Bloch gain profiles 
for lubt — 8 and different values of the magnetic field uIcT =: 
0, 2.5, 3.5. The gain profiles were calculated using Eqs. ([Sj- 



tion of probe field is not possible. By increasing the elec- 
tric field Edc we obtain a completely different situation, 
where high-frequency gain can appear due to nonlinear 
cyclotron oscillations [Fig. [T](d), (e)]. The characteris- 
tic shape of the gain profile turns out to be an inverse 
of the usual dispersive Bloch gain profile [cf. Fig. [T] (d) 
and Fig. [2] (b)]. In this gain profile the frequency of the 
nonlinear cyclotron oscillations Q determines the posi- 
tion of a resonant crossover from loss to gain. When we 
are approaching the separatrix lob = 2wc, we see that 
replicas of this inverted dispersive gain profile appear at 
odd harmonics of Q corresponding to the anharmonicity 
of the ballistic oscillations [Fig. [T] (e)]. With a further 
increase of E^c we arrive to Bloch-like regime ub > Zwc- 
Here we always have a dispersive Bloch gain profile with 
a crossover frequency [Fig. [1] (f )] . Replicas of this gain 
profile can now appear at all harmonics of f2. 

Our results show that both the cyclotron gain and the 
Bloch-like gain have such attractive properties as tunabil- 
ity and extremely large magnitudes of the gain, which can 
be utilized in an operation of THz oscillators and ampli- 
fiers. The gain profiles in both regimes of motion are 
tunable by changing magnetic and electric fields. In the 
Bloch-like regime, the gain profile can be controlled more 
easily with a variation of the electric field, because fl de- 
pends strongly on [see Fig. [T](a)]. In this regime is 
sensitive to the magnetic field only near the separatrix. 
On the other hand, since the cyclotron gain exists only 
near the separatrix, it is tunable by simultaneous varia- 
tion of the electric and magnetic fields as shown in Fig. [2] 
(a). We see that large cyclotron gain can be obtained at 
frequencies lut > f if the cyclotron frequency is some- 
what larger than the scattering rate oJcT > 2 [Fig. [2] (a)]. 
In typical SLs r « 200 fs so that ujt = 1 corresponds to 
the frequency [.i;/27r — 0.8 THz. Therefore the cyclotron 
gain, similarly as the Bloch gain, is ideal for amplifiers 
and oscillators operating at THz frequencies. 

We turn to the important issues of the magnitudes and 
the origin of the high-frequency gain. Fig. [2] (b) shows 
the Bloch gain profiles for fixed ub and different values 
of LUc- We see that by increasing the magnetic field the 



gain increases roughly by an order of magnitude in com- 
parison with the usual Bloch gain at B = 0. Similarly, 
the THz cyclotron gain near the separatrix [Fig. [T] (e)] is 
also an order of magnitude larger than the usual Bloch 
gain. For typical SL parameters d = 6 nm, A = 60 
meV and m — 0.067me, we obtain that uJcT = 1 and 
ujbt = 1 correspond to magnetic field i? = 2 T and 
electric field Edc = 5.5 kV/cm, respectively. Thus, al- 
ready rather weak electric and magnetic fields can pro- 
vide strong THz gain [Fig. [2] (a)] . The magnitude of the 
gain a in units cm~^ is related to the dynamical con- 
ductivity as a = ao(cr,/(To), where at low temperatures 
ao = Ne'^T/cy/eeorrix- For moderate doping N = 10^^ 
cm~'^ and relative permittivity e = 13 (GaAs), we ob- 
tain ao ~ 833.5 cm~^. Using this value of ag, it is easy 
to estimate all magnitudes of the gain in Figs. [T] and [51 
In the vicinity of the separatrix THz gain can have un- 
precedented values a ~ 500 cm^^. 

Although the response to the ac field results from quite 
complicated nonlinear electron dynamics, we can iden- 
tify an important role of the separatrix in the origin of 
this large THz gain. In the Bloch-like regime electrons 
periodically visit the upper part of the miniband and un- 
dergo Bragg reflections at the boundary of the Brillouin 
zone. In contrast there are no Bragg reflections in the 
cyclotron-like regime. However, a strong enough electric 
field Edc periodically brings the carriers to the upper part 
of the miniband where their effective masses are negative 
[11, Hi- We attribute the cyclotron gain to these non- 
classical cyclotron oscillations in presence of weak dissi- 
pation (scattering). As the electrons performing Bloch- 
and cyclotron-like oscillations are spending more time in 
the upper part of the miniband the gain profiles become 
more pronounced. In the vicinity of the separatrix, the 
electron trajectories stick a long time near the Brillouin 
zone boundary (the hyperbolic point of the pendulum) 
resulting in a significant enhancement of the THz gain 
[Figs. [11(c), [21(b)]. 

Since in the cyclotron-like regime the SL operates in 
conditions of PDC, it is inherently stable with respect 
to the space-charge fluctuations. In the gain profiles this 
can be seen as positive dynamical conductivity at low fre- 
quencies [Fig. m (c)-(e), Fig. [2l (a)]. On the other hand, 
in the Bloch-like regime the gain always extends to the 
low frequencies [Fig. [H (f). Fig. [21 (b)] due to the opera- 
tion in the NDC state [Fig. [1] (b)]. This means that in 
the Bloch-like regime the problem of electric instability 
exists similarly as in the case of usual Bloch gain in SLs. 
Several different approaches have been suggested in order 
to stabilize the Bloch gain, including super-superlattice 
quasistatic modulation of the bias |2S 
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structures ^ ^_ ^ ^ ^ 

and 2D shunted surface SLs [13] ■ Here we demonstrate 
that the electric instability can be circumvented by intro- 
ducing an additional magnetic field component in the SL 
direction, so that the magnetic field becomes tilted with 
respect to the SL axis B = {B^, 0, B^)- In this case, elec- 
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FIG. 3: (color online), (a) VI characteristic in tilted magnetic 
field, calculated using Eqs. (|3}-([5]), for UcT — 2 and uic±'r ~ 5 
(black). The thin blue line shows the Esaki-Tsu VI charac- 
teristic [Eq. ((TJ]. (b) Gain profile for the same magnetic field 
and dc bias ujbt — 4.75 (black). The thin blue line indicates 
the usual Bloch gain at B = [Eqs. (|6} and (O]. 



irons perform cyclotron oscillations in the plane perpen- 
dicular to the SL axis with the frequency Wc± = eBx/m. 
These in-plane cyclotron oscillations are coupled to the 
Bloch oscillations via the perpendicular magnetic field 
component Uc- If lob and lUc± are commensurate (Stark- 
cyclotron resonance), the coupling results in delocaliza- 
tion of the electrons, which reveals itself as additional 
resonant structures in the VI characteristics [1]. Fig. [3] 
(a) demonstrates such additional structure for a partic- 
ular magnetic field. The enhancement of the current at 
the Stark-cyclotron resonance lob = Wc_l resembles 
the resonant structures in VI characteristic induced by an 
auxiliary THz field j^. By choosing the working point 
at the PDC part of the peak, we find, using Eqs. ([H])-®, 
that the gain profile close to the Bloch frequency w « 
has the usual shape of the dispersive Bloch gain, whereas 
the dynamical conductivity at low frequencies is now pos- 
itive [Fig.jSjb)]. Interestingly, this stable gain profile can 
be well approximated by a variant of the Tucker formula 
(l6|), in which jdcieEdcd) is now the current density cal- 
culated in the presence of the magnetic field. We should 
distinguish our results from the recent proposal of sub- 
THz generation in SLs in tilted magnetic field 2j| . There 
the moving charge domains are responsible for the cur- 
rent oscillations, whereas in our case the THz gain occurs 
in the absence of electric domains. 

In summary, we found that the magnetic field signif- 
icantly alters the shapes of gain profiles and the mag- 
nitude of THz gain in SLs. We described a novel type 
of large and tunable THz gain caused by nonlinear cy- 
clotron oscillations in the crossed electric and magnetic 
fields. Since the operation point can be chosen at the 
PDC part of the VI characteristic the old problem of 
space-charge instability, which is typical for the Bloch 
gain in SL, does not exist here. We also predicted an 
enhancement of the Bloch gain due to nonlinear charac- 
ter of Bloch oscillations in the presence of a perpendicu- 
lar magnetic field. Finally, we demonstrated that in the 
tilted magnetic field configuration the usual Bloch gain 
can be realized in conditions of PDC. 



We conclude with two remarks. First, by numerically 
solving the Boltzmann equation ^ , we observed that the 
cyclotron gain decreases rapidly with increasing temper- 
ature. It is typically smaller than the Bloch gain [3, 
already at 100 K. The stable THz gain in tilted mag- 
netic field, by contrast, declines slowly and has similar 
magnitudes as the usual Bloch gain also at room tem- 
perature. Secondly, the large-signal response near the 
separatrix can be quite different from the case of linear 
response. In the absence of scattering the electron dy- 
namics can even become chaotic. Nevertheless, we found 
that the cyclotron gain in conditions of PDC can still ex- 
ist also for reasonably large amplitudes of the probe field 
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